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Abstract 



We show for a broad class of counting problems, correlation decay (strong spatial mixing) 
implies FPTAS on planar graphs. The framework for the counting problems considered by us is 
the Holant problems with arbitrary constant-size domain and symmetric constraint functions. 
We define a notion of regularity on the constraint functions, which covers a wide range of natural 
and important counting problems, including all multi-state spin systems, counting graph homo- 
morphisms, counting weighted matchings or perfect matchings, the subgraphs world problem 
transformed from the ferromagnetic Ising model, and all counting CSPs and Holant problems 
with symmetric constraint functions of constant arity. 

The core of our algorithm is a fixed-parameter tractable algorithm which computes the 
exact values of the Holant problems with regular constraint functions on graphs of bounded 
treewidth. By utilizing the locally tree-like property of apex-minor-frec families of graphs, the 
parameterized exact algorithm implies an FPTAS for the Holant problem on these graph families 
whenever the Gibbs measure defined by the problem exhibits strong spatial mixing. We further 
extend the recursive coupling technique to Holant problems and establish strong spatial mixing 
for the ferromagnetic Potts model and the subgraphs world problem. As consequences, we have 
new deterministic approximation algorithms on planar graphs and all apex-minor-free graphs 
for several counting problems. 

1 Introduction 

In study of counting algorithms, many counting problems can be formulated as computing the 
partition function: 



where <&e '■ [q] 2 —> C and <£y : [q] — > C are symmetric functions. This model is called spin system 
in Statistical Physics. It has vertices as variables and edges as constraints, and the partition 
function returns the total weight of all configurations. Many natural combinatorial problems such 
as counting independent sets, q-colorings, or graph homomorphisms can be expressed in this way. 

We consider a framework that encompasses a much broader class of counting problems, namely, 
the Holant problems. 

* Supported by the National Science Foundation of China under Grant No. 61003023 and No. 61021062. 
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An instance of a Holant problem is an Q = (G(V,E), {fv}vev), where G is a graph, and each 
f v is a function that maps tuples in [<7] de s( u ) to function values. The Holant of is defined as 

hoi(o)= ]T Y[f v {*\ E{ v)), 

ae[q] E v£V 

where f v {a \e(v)) evaluates f v on the restriction of a on incident edges E(v) of vertex v. The Holant 
problem Holant (Q, T) specified by a graph family Q and a function family J 7 , is the problem of 
computing hol(O) for all valid instances defined by graphs from Q and functions from T . 

The term Holant is coined by Valiant in |56j in studying of holographic algorithms. The formal 
framework of Holant problems is proposed in [TTJ by Cai, Lu and Xia. 

The Holant framework is extremely expressive. Using the bipartite incidence graph to represent 
the participants of variables in constraints and choosing appropriate functions at vertices on both 
sides, computing the partition functions of spin systems and more generally counting CSPs can all 
be represented as special classes of Holant problems. 

An algorithmic significance of Holant problems is that they are outcomes of holographic trans- 
formations. The holographic algorithms proposed by Valiant |55p56j compute exact solutions to 
the counting problems on planar graphs by transforming to problems solvable by the FKT algo- 
rithm [25U38tl5^] for counting planar perfect matchings. In the realm of approximate counting, 
perhaps the most successful (implicit) using of holographic transformation and Holant problem is 
the FPRAS for ferromagnetic Ising model given by Jerrum and Sinclair in their seminal work [36J. 
The transformation in [36] from the spins world to the subgraphs world is indeed a holographic 
transformation, and the resulting subgraphs world problem is a Holant problemQ In these exam- 
ples, the original counting problem is transformed to a Holant problem which has efficient exact or 
approximate algorithms. Therefore, the following problem is fundamental to the study of counting 
algorithms: 

Characterize the tractability of exact computation and approximation of Holant (Q, J-) in terms 
of graph family Q and function family T . 

Exact computation. The exact computation of Holant problems has been well studied on gen- 
eral graphs [7ll8lfi m[T3Tfl5ll43"] and planar graphs [9l[l2j[56], sometimes in form of dichotomy the- 
orems, which states that every problem in the considered framework is either ^P-hard or having 
polynomial-time algorithm. Very recently, a dichotomy theorem [7] is proved for Holant problems 
with complex- valued functions on general graphs, concluding a long series of dichotomies on Holant 
problems. All these results consider Holant problems with boolean domain (q = 2). Meanwhile, 
some special classes of Holant problems are more thoroughly understood, such as counting graph ho- 
momorphisms or counting CSP. For these specialized frameworks, dichotomy theorems were proved 
in a very general setting with complex- valued functions on general-sized domains (6j[16] . See |18| 
for a good survey on these subjects. 

Speaking very vaguely, the dichotomy theorems tell us that except for some rare cases almost 
all Holant problems are hard. Then a problem of algorithmic significance is to establish tractable 
results for Holant (Q, F) on more refined graph families Q, e.g. graphs with fixed parameters or 
forbidden minors, especially for general domain size q > 2. 

lr This actually happened more than a decade earlier than the concepts of holographic algorithm and Holant 
problem formally defined. 
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Approximation. We focus on deterministic approximate counting algorithms, specifically, the 
deterministic fully polynomial time approximation scheme (FPTAS). A central topic in this di- 
rection is the relation between correlation decay (strong spatial mixing) and approximability of 
counting. 

Correlation decay is a property about the marginal distribution, which is computationally equiv- 
alent to counting by the famous self-reduction of Jerrum- Valiant- Vazirani [37]. The correlation 
decay property says that faraway vertices have little influence on the marginal distribution of local 
states, thus marginal probabilities should be well-approximated by local information only. However, 
as noted in [5 , 33] , this sufficiency of local information does not immediately yield efficient local 
computation. Two tools are invented to bridge this gap: the self-avoiding-walk tree (SAW-tree) of 
Weitz |57j and the computation tree of Gamarnik and Katz [33| . Both transform the original graph 
to a tree structure in which the marginal probabilities can be efficiently computed by recursions. 
With the SAW-tree the implication from strong spatial mixing to FPTAS is proved for 2-state spin 
systems |57] . which becomes a foundation for several important algorithmic results |41U42[ |4"8 ] I50]. 
It is also proved in a long series of beautiful work [23][3lJ[32j|47][!2l[53] that for the same class of 
counting problems lack of correlation decay implies inapproximability. 

The relation between correlation decay and approximability for broader classes of counting 
problems is widely open, because of following technical challenges: 

• It is known [51] that for domain size q > 2 tree may not always represent the extremal case 
for correlation decay. Thus in order to use correlation decay to support approximate counting 
for those problems, the local computation has to be done on structures other than trees. 

• Even on trees, the current recursion-based computation critically relies on the simplicity of 
constraint functions, as in the cases of spin systems and matchings. For general Holant 
problems, even on trees and when q = 2, it is not known whether simple recursion exists. 

1.1 Our results 

We make progress on both exact and approximate computation of Holant problems by establishing 
connections between them. 

We characterize a broad class of Holant problems whose exact computation is tractable on 
tree- like graphs and FPTAS is implied by strong spatial mixing on planar graphs. These Holant 
problems are characterized by a notion of regularity introduced by us on the constraint functions. 
Intuitively, being regular as a function means that the entropy of any input or partial input is 
constant. This covers a large family of important counting problems, including all spin systems, 
graph homomorphisms, counting CSP with symmetric constraints of bounded arity, matchings, 
perfect matchings, the subgraphs world problem in [36], etc. 

For this broad class of counting problems, we give a fixed-parameter tractable algorithm which 
computes the exact value of counting in time 2°^ ■ poly(ra) on graphs of size n and treewidth k. 
Based on this parameterized algorithm, strong spatial mixing implies FPTAS on apex-minor-free 
graphs, which include planar graphs as special case. 

We also apply the recursive coupling technique of Goldberg et al. [35] to analyze the strong 
spatial mixing for Holant problems. As examples, we have deterministic FPTAS on planar graphs, 
and more generally on all apex-minor- free graphs for the following counting problems: 
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• Counting g-colorings on triangle-free planar graphs of maximum degree A when q > aA — 7 
where a ~ 1.76322 and 7 ~ 0.47031. This is just directly applying |35j . 

• The subgraphs world with parameter fi, A < 1 on planar graphs of maximum degree A when 
A < ! and as a consequence the ferromagnetic Ising moderj with inverse temperature 

(3 and external field B when A < ^b^+I^Jb ■ 

• Ferromagnetic q-state Potts model of inverse temperature (3 on planar graphs of maximum 

j / q — 2 \ 

degree A when /5 < — a+i 1 , which vastly improves the mixing condition in [33] for FPTAS 
on general graphs and is close to the j3 = O(^) bound conjectured in [33]. 



Technical contributions. Our parameterized algorithm does not directly use the tree decompo- 
sition. Instead, we define a new decomposition called the separator decomposition, which recursively 
separates the graph by small graph separators into components of limit-sized boundaries. This is 
quite different from the known treewidth-based approaches for spin systems, e.g. the junction tree 
algorithm; and this new construction more closely aligns with the conditional independence prop- 
erty: conditioning on any fixed assignment on a separator, the states of separated vertices are 
independent. The construction of separator decomposition makes explicit connections between the 
separable structure of tree-like graphs and the conditionally independent nature of counting prob- 
lems defined by local constraints. As a result, our algorithm can deal with much broader class of 
counting problems other than just spin systems. 

Unlike previous approximation algorithms via correlation decay, where the decay is verified on 
a tree of size exponential in the size of original graph, our FPTAS only relies on the correlation 
decay on the original graph. Thus we can directly apply those "decay-only" results such as [35] to 
get FPTAS. Since we do not explode the size of the graph, the FPTAS can even be supported by 
single-site correlation decays. 



1.2 Related work 

The use of correlation decay technique for designing FPTAS for counting problems was initiated 
in [H[57] and has been successfully applied to many problems [3[3T1[33], especially for computing 
the partition function of Ising model [41U421 I50] . The technique of recursive coupling has been used 
to prove the property of correlation decay [MIESEIIIIH] • 

The locally tree-like property of planar graphs and apex-minor-free graphs provides structure 
information to develop both exact and approximation algorithms on decision and optimization 
problems, some examples include [3"ll2T | l24" l [2"T] ■ 

A framework for parameterized complexity of counting problems was proposed in [2l l28ll4^] . 
The parameterized complexity of computing partition functions has been studied via probabilistic 
inference in graphical model [17] . Some logical approaches have also been extended to counting 
problems on structures with small treewidth or local treewidth [TJdHlEO]. 

2 Due to the classic results of 1960s in Statistical Physics [26139] and the recent theory of holographic algorithms [9], 
the planar Ising model with zero field B = is solvable exactly in polynomial time. Here we consider Ising models 
with general field B, which is #P-hard even on planar graphs proved implicitly in [40] . 
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2 Models and statement of results 



2.1 Holant problems 

Let [q] = {0, 1, . . . , q — 1} be a domain of size q, where q > 2 is an finite integer. Let / : [q] — > F 
be a d-ary function where F is a field. In this paper, we consider either F = C the complexes or 
F = R + the nonnegative reals. To avoid issues of computation model, we assume all number are 
algebraic. We allow the function arity d to be 0. When d = 0, the only member of [q]° is the 
empty tuple £, and a 0-ary function / maps £ to a function value. We call such function / a trivial 
function. 

A d-ary function is symmetric if f(x±, . . . ,Xd) = f(x p m, ■ ■ ■ > x P (d)) f° r an Y permutation p of 
{1, 2, . . . , d}. When q = 2, functions have boolean domain and a d-ary symmetric function / can be 
denoted by [/o, /i, ■ ■ ■ , /d] where specifies the function value for the input tuple with Hamming 
weight k. For example the Equality function is denoted as [1, 0, ... , 0, 1]. 

Let $e '■ [<l] 2 C and <3?y : [q] — > C be two symmetric functions. The partition function of an 
undirected graph G(V, E) is defined as 

Z ^ = E II *e(<t(u),o-(v))H$ v (o-(v)). 

ae[q] v {u,v}eE vev 

This is called a q-state spin system. Let Q = (G(V, E), {f v }vev) be an instance, where each f v , 
called a constraint function or a signature, is a d-axy symmetric function with d = deg(v). We 
define the Holant of f2 as 

hoi(n)= ]T n^^l^))' 

where f v {a \e( v )) evaluates f v on the restriction of a on incident edge of v. 

Let Q be a family of graphs and J 7 be a family of functions. A Holant problem Holant(t?, J-) is 
a computation problem that given as input an instance £1 = (G(V,E), {f v }v&v) where G € Q and 
all f v are from J 7 , compute hol(Jl). 

A symmetric function can be represented by a vector enumerating the function values for all 
inputs (up to symmetry). The number of symmetry classes of a & [q] d equals the number of weak 
(/-composition of integer d, which is ( d+< ^ 1 ) ■ Thus symmetric d-ary functions can be represented 
by vectors of length polynomial in d. 

Spin systems can be represented as special class of Holant problems. For a graph G(V,E), let 
Tq denote the incidence graph of G, i.e. Tq = (Vi, V2, E') is a bipartite graph with V\ = V, V2 = E 
and (v, e) € E' if and only if edge e is indecent to vertex v in G. For a spin system defined by 
functions &e and &v on a graph G, we can transform it to a Holant instance 17 — {ILq, {fv\v£VuE) j 
where f v = <&e for right vertices v £ E and for left vertices v G V, f v is the generalized Equality 
function defined as f(x±, . . . , xj) = &v( x i) if x i = " " " = x d an d f(xi, ■ ■ ■ , x^) = if otherwise. It 
is easy to check that hol(ri) = Z(G). 

2.2 Regular functions 

Our characterization of Holant problems relies on a "pinning" operation on symmetric functions. 
The pinning operation on a function defines a new function with smaller arity by by fixing (pinning) 
the values of some of the variables. 
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Definition 2.1 (pinning). Let f : [q] d — > F be a d-ary symmetric function. Let < k < d and 
t G [q] k - We define that PlN T (/) = g where g : [q] d ~ k — >■ F is a (d — k)-ary symmetric function 
such that 

Va G g{a) = f(a(l), ...,a(d- k), r(l), . . . , r(fc)). 

Specifically, when k = i/ie resulting function g = f ; and when k = d, the resulting function g is a 
trivial function f(o~). 

Note that since / is symmetric, the positions of r(l), . . . ,r(k) in / does not matter, and the 
pinning of a symmetric function is still symmetric. 

To exemplify the effect of pinning, consider the case when q = 2 and a function / is represented in 
form [f , fi,...,fd]. For a a G [2} k that a has I many Is, we have PlN ff (/) = [f t , f e+1 ,. . . , fd-(k-i)]- 
That is, the Pin ct (/) for a a G [2] k returns a "sliding windows" of length d — k in [/o, fi, ■ ■ ■ , fd] 
whose starting position is determined by the number of Is in a. 

A notion of regularity of symmetric functions can be defined by limiting the outcomes of pinning. 

Definition 2.2 (constant regularity). A symmetric function f : [q] d — > F is called C-regular if for 
all < k < d, it holds that 



{pin t (/)|Vtg[#} 



< C. 



A family T of symmetric functions is called regular if there exists a finite constant C > such that 
every f G T is C-regular. 

We then give some sufficient conditions for regular symmetric functions. 

Proposition 2.3. Let f : [q] d — >■ F be a symmetric function. For a G [q] and i G [q], let 
n>i(o~) = |{1 < j < d | o~(j) = i}\ be the number of i- entries in a. 

• (bounded arity) f is ( d ~^ ( ^ 1 ) -regular. 

• (cyclic) If there is a c > such that f(o~) depends only on {n\{o) mod c, . . . ,n q (a) mod c) 
then f is c 9_1 -regular. 



• (constant exceptions) If there is a C-regular g : [q] d —> F and a c > such that f and g differ 
only at those a G [q] d that nj(cr) > d — c for some i G [q], then f is [c + q ■ i^^^ 1 ^ -regular. 

Therefore all constant-ary symmetric functions, Equality and generalized Equality, and the 
Not- All-Equal are all regular. This covers all g-state spin systems. 

For boolean domain, a function [/o, fi, . . . , fa] is regular either if it is cyclic, i.e. fk = mo d c 
for some constant c, or if it becomes cyclic after removing constant many exceptions /o,/i, • • • ,/c 
and fd-c, ■ ■ ■ 5 fd from both ends. This covers (weighted) matchings, perfect matchings, and the 
subgraphs world transformed from the Ising model [86] , which is a Holant problem defined by 
constraint functions in the form [1, fi, 1, fi, . . .] and [1, 0, A]. 



2.3 Correlation decay 

Let = (G(V, E), {fv}vev) be a Holant instance where each f v : [q\ de s( v ) — >. U+ i s a symmetric 
function with nonnegative real function values. Let a G [q] E be a configuration and w(a) = 
Yivev fv i a \e(v)) be the weight of configuration a. 
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A configuration a € [q] E is feasible if w(a) > 0. And for a configuration ta € [q]^ on a subset 
A C E of edges, we say that ta is feasible if there is a feasible a G [gj^ agreeing with ta on A. 

The Gibbs measure is a probability distribution over all configurations, defined as /i(<t) = g^fe^ ■ 
To make the Gibbs measure well-defined, we require that each /„ has nonnegative values and the 
Holant problem is feasible, i.e. there exists a feasible configuration. 

For a feasible a\ £ [gf on A C E, we use [i a e K to denote the marginal distribution at e 
conditioning on the configuration of A being fixed as cr\. 

Definition 2.4 (Strong Spatial Mixing). A Holant problem Holant((/, J 7 ) has strong spatial mixing 
(SSM) if for any instance Q, = (G(V,E),{f v } ve v), o- n V e £ E, A C E and any two feasible 
configurations oa,ta € it holds that 

\K A - mI a IItv < Poiy(l^l) • exp(-n(dist( e> A))), 

where A C A is the subset on which o~\ and t\ differ, dist(e, A) is the shortest distance from edge 
e to any edges in A, and \\ ■ ||tv denotes the total variation distance. 

The strong spatial mixing defined for spin systems [57] is covered as special case. 

2.4 Tractable search 

In order to apply the self-reduction technique of Jerrum- Valiant- Vazirani [37] for approximate 
counting, we also require that the following search problem is tractable: 
Input: a Holant instance Q, = (G(V, E), {f v } v ev)> and a configuration a\ € [q] A on A C E; 
Output: a feasible r £ [q] E agreeing with oa on A, or determines no such r exists. 

We call such property the tractable search for Holant (Q, F). We remark that this is a very 
natural assumption for approximate counting: for all known examples of approximate counting im- 
plied by mixing, the above search problem is easy or even trivial. The tractable search requirement 
of the general Holant framework is an analog to the specific q > A + 1 requirement for counting 
g-coloring. 

The tractable search is related to the Holant c framework. For i £ [q], let Aj denote the unary 
function which maps i to 1 and all other j € [q] to 0. 

Definition 2.5. Holant c (£, F) = Holant (Q, T U {Aj | i € [q]}). 

The Holant c problem has significance in complexity of counting [El|13j . Assuming the tractable 
search for Holant(£/, J 7 ) is equivalent to assuming the polynomial-time decision oracle (for existence 
of feasible configuration) for Holant c (<?, J 7 ). 

2.5 Local treewidth and planarity 

Our algorithm relies on the treewidth of graph and family of graphs with forbidden graph minors. 
We will not formally define these concepts excepting saying that the treewidth measures how similar 
a graph is to a tree. The formal definitions can be found in standard textbooks, e.g. [22] . 

Graphs of bounded local treewidth are precisely the family of apex-minor-free graphs, where an 
apex graph has a vertex whose removal leaves a planar graph. In particular, K$ and .^3,3 are apex 
graphs, therefore apex-minor- free graphs include planar graphs as a special case. 

Let tw(G) denote the treewidth of graph G. 
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Theorem 2.6 ( [20^124] ) . Let Q be an apex-minor-free family of graphs. For any G(V,E) G Q and 
v G V , let N r (y) be the subgraph of G induced by vertices whose distance to v is at most r. Then 
tw(N r (v)) < f(r) for some linear function f. 

The following easy proposition states that representing spin systems as Holant problem on the 
incidence graphs does not violate the graph structure. 

Proposition 2.7. Let G be a graph and Xq be the bipartite incidence graph of G. Then tw(G) = 
tw(Xc) and Tq is apex-minor-free if G is apex-minor- free. 

2.6 Main results 

Our main results can be summarized by the following two theorems. 

Theorem 2.8. There is an algorithm for Holant (Q, F) with regular symmetric T , whose running 
time is 2°^ ■ poly(n) for any G G Q of n vertices and treewidth k. 

Theorem 2.9. Assuming the tractable search, for Holant problem Holant (Q, J 7 ) with apex-minor- 
free Q and regular symmetric nonnegative T , SSM implies FPTAS. 

3 Structure of regular functions 

In this section we introduce a new concept which characterizes the structure of outcomes of pinning 
as well as plays a key role in efficient computation of Holant problems. 

Note that different a, r G [q] k (up to symmetry) may yield the same function after pinning / 
with them. We classify members of [q] k into equivalence classes according to their effects of pinning 
on / by introducing the following concept of peers. 

Definition 3.1 (peers). Let f : [q] d — > ¥ be a d-ary symmetric function. Let < k < d and 
t G [q] k ■ We define that Peer t (/) = g where g : [q] k — > {0,1} is a boolean symmetric function 
such that 

w cMfc , s fl ifPlN CT (/)=PlN T (/), 
I otherwise. 

We may also interpret Peer t (/) as a set and write a G Peer t (/) if Peer t (/) (a) = 1. 

The condition Pin ct (/) = Pin t (/) defines an equivalence relation between a and r (having the 
same effect of pinning on /). Then Peer t (/) is the indicator function of the equivalence class 
{a G [q] k | Pin ct (/) = Pin t (/)}. So we have the following easy but useful proposition. 

Proposition 3.2. a G Peer t (/) if and only if Peer ct (/) = Peer t (/). 

The peer images of an input uniquely determines the function value, specifically: 

Lemma 3.3. Let f : [q] d —> ¥ be a symmetric function. Let r > 1, d\ + c?2 + ■ ■ ■ + d r = d, 
and ai,Ti G [q] di for i = 1,2, ...,r. If Peer^. (/) = Peer t . (/) for all i = 1,2, ...,r, then 
/(tit 2 • ■ ■ T r ) = f(aia 2 ■ ■ ■ cr r ). 
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Proof. Due to Proposition 13.21 r, G Peer^. (/) for all i = 1,2, ... ,r. We prove by induction on r. 
When r = 1, for all d > and o"i € [g] rf , Pin (/) is a trivial function f(o~i) (a function value) 
and Peer (7i (/) is the equivalent class of all t\ G [q] d which have the same Pin (/) as Pin (/), 
i.e. f{n) = /((Ti). 

Assume the statement holds for all smaller r and all d. Let <7j, Tj € i = 1,2, ... ,r satisfy 
that n G Peer,,. (/) for all i = 1, 2, . . . , r. Since r r € Peer ct? . (/), we have PiN Tr (/) = PiN CTr (/). 
Denote that g = PlN Tr (/) = PlN CTr (/). By definition of pinning, it holds that f{p\Ui ■ ■ • oy) = 
g{o~\o~2 ■ ■ ■ o>- 1) an d /(T1T2 • • • T r ) = gij\T2 ■ ■ ■ T r -\). Note that g satisfies the induction hypothesis 
for r — 1, which means that g{r\T<i ■ ■ - r r _i) = g(o\o~i ■ ■ ■ 0>_i) as long as Tj G Peer (/) for 
all % = 1,2, . . . ,r — 1. Therefore, for <7j,Tj € [<?] d V = 1,2, ...,r that Tj G Peer ct . (/) for all 
i = 1,2, ... ,r, we have /(rir 2 • • • r r ) = g(nT 2 ■ ■ ■ r r _i) = 5(cricr 2 • • • oy-i) = /(cri£T 2 • • • oy). □ 

Note that the outcome of Peer t (/) is still a symmetric function, so we can apply pinning and 
peering operations on it. We can define the peering closure which contains all possible outcomes of 
recursively applying peer operations on a function /. 

Definition 3.4 (peering closure). Let f : [q] d — > ¥ be a d-ary symmetric function. Let < k r < 
k r -\ < ■ • • < k± < d and Tj G [q] ki , 1 < i < r. We denote that 

PEEB n,7S,...,r P (/) = PEE Rr r (PEER 71j7a> ... iTr _ 1 (/)) . 

The peering closure of f, denoted by Peer* (/), is defined as 
Peer* (/) = {peer T1iT2| ... i7V (/) | n G [qp for 1 < i < r, < k r < k r -i < ■ ■ ■ < h < d, r > l} . 

Note that Peer t (/) is a boolean functions no matter what the range of / is. A boolean function 
g : [q] d — > {0, 1} can be seen equivalently as a set {a G [q] d \ g(o~) = 1}. For two boolean functions 
g and h defined on the same domain \q] d , we define the operations on boolean functions gUh, gC\h, 
and g C h according to the operations on their set representations. 

Some useful properties of peering are better presented in this set language: 

Lemma 3.5. Let f : [q] d — > F be a symmetric function, and g, h : [q] d — > {0, 1} be boolean symmetric 
functions. Let a G [q] and r G [q] k for arbitrary < k < t < d. We have 

1. Peer t (/) C Peer t (Peer, (/)); 

2. (Peer t (g) n Peer t (h)) C Peer t (gUh). 

Proof. Both statements can be proved by directly expanding the definition of Peer (•). 

1. For any ir G [q] h , suppose that ir G Peer t (/), which means Pin,,. (/) = PlN T (/). Then 
for any x G [qf~ k , we have PlN^ (/) = PlN XT (/), thus PlN X7r (/) = Pin, (/) if and only 
if PlN xr (/) = PlN CT (/), which is equivalent to that for any x G [qf~ k , xir G Peer ct (/) if 
and only if xt G Peer ct (/). This implies that Pin,,. (Peer ct (/)) = PiN r (Peer ct (/)), which 
implies n G Peer t (Peer ct (/)). Therefore, Peer t (/) C PEER r (Peer ct (/)). 

2. For any tt G [q] k , suppose that tt G Peer t (g) n Peer t (h), which implies that Pin t (g) = 
PlN^. (g) and PlN r (h) = PlN^ (h). Then for any x G [q] d ~ k , it holds that xt G g if and only if 
xtt G g and xt G h if and only if xtt G h, thus xt G gUh if and only if xtt G gUh, which means 
PiN r (gUh) = PiN n {g U h), thus it G Peer t (gUh). Therefore, (Peer t (g) n Peer t (h)) C 
Peer t (gUh). 
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□ 

We then characterize all k-ary functions in peering closure by unions of boolean functions. 

Lemma 3.6. Let f : [q] d — > F be a symmetric function. For any < k < d, every k-ary function 
g G Peer* (/) can be represented as 

t 

g = |J PEER CTi (/) 
i=i 

for some ai,... ,a t G [q] k . 

Proof. Recall that every fc-ary function in Peer* (/) is in the following form: 

PEER Tl , T2 ,..., Tr (/) = PEER Tj . (PEER Tr l (• • • (PEER Tl (/)) • • 

for some r > 1, k = k r < k T ~\ < • • • < k\ < d, and T{ G [q] ki , 1 < i < r. 

We then prove by induction on r that Peer ti T2 Tr (f) = Peer (Ti (/) U • • • U PEER CTt (/) for 
some <7i, . . . ,<Tj G [q] k where k is the arity of r r . When r = 1, this is trivially true. Assume the 
statement holds for all smaller r. Then Peer ti T2j Tr (/) = PEER Tr ^PEER T1 )7V _ x And due 

to the induction hypothesis, there exist a\, . . . ,o~t G [q]^- 1 such that 

t 

Peer^...^ (/) = Jpeer CTi (/). 

i=l 

Due to Lemma 13.51 we have 

t / t \ 



p| PEER Tr (PEER CTi (/)) C PEER Tr (J PEER CT . (/) ; and 
i=l \i=l J 

PEER Tr (/) C PEER Tr (Peer ct (/)) for any a G for]*- 1 . 

Combining these, we have 

PEER Tr (/) C PEER Tr ( (J PEER^ (/)) = PEER Tr (PEER^...^ (/)) = PEER ri)r2) ... )rr (/) . 



\i=l 



Note that this already implies the lemma, i.e. there exist t > 1 and ai,...,o~t G [q] k where k 
is the arity of r r such that Peer ti T2j Tr (/) = (J i=1 f >EER o- i (/)• To see *hi s implication, by 
contradiction we assume that the statement is false. Then there must exist a\,a 2 G [q] k such 
that o-i,a 2 G Peer cti (/) but cii G Peer tiT2j Tr (/) and cr 2 Peer ti T2j Tr (/). Recall that 

Peer^ (/) are equivalent classes. Then a± G Peer ti T2 )7V (/) = Peer Ti . ^Peer t1i Tr _ 1 (/)^ 
implies that Peer^ ^Peer^ r „ )Tr _ 1 = Peer t ^ ^Peer^ r .. )Tr _ 1 (/))• On the other hand, we 
have shown that Peer cti (/) C Peer (Ti ^Peer ti Tr _ ± (f)j . However, due to the assumption 
we have a 2 PEER n>T2) Tr (/) = Peer^ (Peer t1j Trl (/)J, which implies that Peer cti (/) % 
Peer^ ^Peer t1j jTr l (/)), a contradiction. 

□ 
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Lemma 3.7. Let f : [q] d — > F be a symmetric function. If f is C -regular, then 

1. for any < k < d, the number of distinct k-ary functions in PEER* (/) is at most 2 C ; 

2. for every g G Peer* (/), g is C -regular. 

Proof. First note that for any < k < d, it holds that 

{PEER r (/) | T € [q] k }\ = |{PlN r (/) | T G [q] k } . 

This is because they both count the number of equivalence classes defined by the equivalence 
relation that a and r are equivalent if and only if PlN CT (/) = PlN T (/). If / is C-regular, then 
|{Pin t (/) | re [q] k } \ <C, so we have |{Peer t (/) |tG [q] k } \ <C. 

Due to Lemma 13.61 every /c-ary function in Peer* (/) can be represented as a (boolean- 
functional) union U* =1 Peer ct . (/) for some ui,... ,o~t G [q] k . The number of different unions is 
obviously bounded by the size of power set of {Peer t (/) | r G [q] k }, which is at most 2 C . There- 
fore the number of distinct fc-ary functions in Peer* (/) is at most 2 C . The first part of the lemma 
is proved. 

Let < k < d and r G [q] k . By definition of Peer (•), for any < I < k, both {Peer ct (/) | a G [q] e ] 
and {Peer^ (Peer t (/)) | a G [q\ £ } are partitions of [q] e . Due to Lemma 1331 for any a € [q] e , we 
have 

Peer,, (/) C Peer ct (Peer t (/)) . 

Thus {Peer^ (/) | a G [q] e ] is a refinement of {Peer ct (Peer t (/)) | a e [q] e }- Therefore for any 
< £ < k, we have that 

PIN CT (Peer t (/)) I a € [q] £ }\ = IjPEER, (Peer t (/)) | a G [q] 1 } 

< IPEER CT (/)|aG[^} 



{PlN CT (/)|aG[^} 



which is bounded by C if / is C-regular. Therefore, Peer t (/) is C-regular if / is C-regular. The 
second part of the lemma is a trivial consequence of this fact. 

□ 



4 The Separator Decomposition 

In this section we introduce a construction called separator decomposition which is essential to the 
efficient computation of Holant problems. We also show that the width of this decomposition is 
asymptotic equivalent to the treewidth of graphs and give an parameterized algorithm to construct 
separator decomposition. 

Let G(V, E) be an undirected graph and H(U, F) be a subgraph of G. For any vertex set it! C V, 
the vertex boundary of R in H, denoted by dnR, is defined as d}jR = {u £ U\R \ 3v G R,uv G F}. 
In particular, when H = G we omit the subscript and denote the vertex boundary by OR. 

Definition 4.1 (Vertex Separator). Let G(V,E) be an undirected graph. A vertex set S C]/ is an 
(X, y)-separator in G, where X, Y C V are two vertex sets, if 



11 



1. {X, Y, S} is a partition ofV; 

2. all u € X are disconnected from all v £7 in G[V \ S]. 

Note that in the above definition we do not require X and Y to be connected or even nonempty. 
Any vertex set is an (X, Y)-separator if one of X, Y is empty. 
We now introduce the separator decomposition. 

Definition 4.2 (Separator Decomposition). Let G(V,E) be an undirected graph. A separator 
decomposition of G is a full binary tree T (a rooted tree in which each node either has two children 
or is a leaf) with each node iofT associated with a pair (Vi, Si) such that Vi,Si C V and satisfy: 

1. V r = V for the root r in T and Vi = Si = for each leaf I in T; 

2. for every non-leaf node i and its two children j, k in T , Si is a (Vj, Vk)-separator in Gi, where 
Gi = G[Vi] is the subgraph of G induced by Vi. 

The width of a separator decomposition is the maximum of \dVi\ and \Si\ over all nodes i in T. 

For a graph G of n vertices, a separator decomposition T must have O(n) nodes, because for 
each node i € T and its children j, k, we have Vj and disjoint and Vj, C V,. 

We then connect the width of separator decomposition to the well-known treewidth of graphs, 
and give an algorithm for constructing separator decomposition with small width if there is such a 
decomposition for the input graph. 

Theorem 4.3. Let sw(G) be the minimum width of all separator decompositions of G, and let 
tw(G) be the treewidth of G. Then 

I ■ tw(G) + — < sw(G) < 3 • tw(G) + 3. 
8 24 

And there is an algorithm which given as input a graph G of n vertices and treewidth k constructs 
a separator decomposition of width at most 3(k + 1) in time 2°^ ■ poly(n). 

It is possible to have an 0(2?^ ■ n) time algorithm for super-linear f(k), which has better 
performance in the sense of parameterized complexity. However, our FPTAS on apex-minor-free 
graphs critically relies on that the time growth in treewidth k is 2°^ k \ 

We will not give the formal definition of treewidth, but instead we will use the following notion 
of balanced separators to characterize treewidth . 

Definition 4.4 (Balanced separator). Let G(V,E) be an undirected graph and W cy '. A vertex 
set S C V is an balanced P^-separator in G if W \ S can be partitioned into X and Y such that 
< \X\, \Y\ < §|W|, and all u € X are disconnected from all v G Y in G[V \ S]. 

A relation between treewidth and balanced separator is stated in next theorem implicit in |49j . 

Theorem 4.5 (Robertson and Seymour |49j). Let G = (V, E) be an undirected graph. Iftw(G) < k, 
then for every W C V of size at least 2k + 3 there exists a balanced W -separator of size at most 
k + 1. Conversely, if for every W cy of size 3k + 1 there exists a balanced W -separator of size at 
most k + l, then tw(G) < 4k + 1. 

We then use balanced separators to characterize the width of separator decompositions. 
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Lemma 4.6. Let G = (V, E) be a graph of n vertices. 



1. If G has a separator decomposition of width s, then for every W C V of size at least 6s there 
is a balanced W -separator of size at most 2s. 

2. If for every W C V of size 6s there is a balanced W -separator of size at most 2s, then G 
contains a separator decomposition T of width at most 6s. And such T can be constructed in 
time 2 ( s ) • po ly( n). 

Proof. We show the first part: existence of a separator decomposition of width s implies that for 
every W C V there is a balanced ^-separator of size at most 2s. 

Let T be a separator decomposition of G of width s, with each node i S T associated with a 
vertex set Vi and a separator Si of G[Vj]. It holds that \dVi\ < s and \Si\ < s for all i £ T. 

Fix an arbitrary W C V with \W\ > 6s. Let i be the node in T of maximum depth satisfying 
\Vi Pi W\ > t}\W\. Such node i always exists and must be a non-leaf since V r fl W = W for the root 
r and V e n W = for every leaf £. Let 5 = 8Vi U It holds that |5| < + \Si\< 2s. We then 
show that S is a balanced ^-separator. 

Let j, k be the children of i in T. Due to the maximality of the depth of node i, it holds that 
\Vj nW\ < \\W\ and \Vk n W\ < ^\W\. We distinguish between the following two cases: 

Case.l: \VjC\W\ < \\W\ and \V k nW\ < ^\W\. Let X = (VjnW)U(V k nW) and Y = W\{XUS). 
We have \X\ < |V,-rW| + |y fe rW| < l\W\ and \Y\ = \W\{XUS)\ < |W|-|^nW| < \\W\ < ||W|. 
Moreover, both X and V are nonempty, since 



It then remains to show that X and Y are separated by S in G. It is easy to verify that X C Vi 
and y n Vi = 0, thus every path X to Y must go through some vertex in dVi C S. 

Case.2: \Vj n W| > i|W| or |V fe n W| > ||W|. Without loss of generality, suppose that 
\Vj HW\> \V k n W\ and |Yj n W\ > \\W\. Let X = Vj n Y = W\(S U X). We have that 
\X\ = \Vj fl W\ < i|W| and |Y| < |W| - |Yj H W| < §|W|. Also X and Y are nonempty since 
\X\ = \Vj r\W\> \\W\ > and |Y| > \W\ - \X\ - \S\ >\\W\- \S\ > 3s - 2s > 0. 

We then show that X and Y are separated by S in G. It holds that X C Vj and Y n Vj = 0, 
thus every path X to Y must go through some vertex in dVj C (dVi U 5j) = S. 

Therefore, in both cases S is a balanced TY-separator. The first part of the lemma is proved. 

We then prove the second part: if for every W C V of size 6s there is a balanced lY-separator 
of size at most 2s, then G contains a separator decomposition of width at most 6s. We first prove 
the following claim. 

Claim. If for every W C V of size 6s there is a balanced W -separator of size at most 2s, then for 
any nonempty ficy with \dR\ < 6s, there is a partition {X, Y, S} of R such that 

1. S is an (X,Y)- separator in G[R] and \S\ < 4s; 



x\ = \(Vjnw)\Si\ > \vnw 



y| = \w\(xus)\ > \w\ - \x 




2. \dX\,\dY\ < 6s. 
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Proof. When \R\ < 4s, the claim holds by taking S = R and X = Y = 0. We then consider only 
the case \R\ > 4s. 

Let W D dR and \W\ = 6s. Let S' be a balanced VK-separator in G of size at most 2s. Then 
W\S' can be partitioned into X\y and Yyy such that Xy/ and Yy/ are disconnected in G[V\S'], and 
< \X W \, \Y W \ < l\W\S'\. Since G[V\S'] is disconnected, we have that V\S' can be partitioned 
into Xy and Yy such that Xw C Xy, Yw C Yy, AV and Yy are disconnected in G[V \ S'], and 

o < \x v nw\,\Y v nw\ < %\w\s'\. 

We define that S = S' C\ R, X = Xy D R and Y = Yy D R. We then verify that they satisfy the 
requirement. 

If X = R (or Y = R) then 5 = 0, and Xy (respectively Yy) contains more than 4s vertices in 
W, contradicting that 5' is a balanced VF-separator. And it follows from that S' separates Xy and 
Yy in G that S is an (X, Y)-separator in G[R\. It trivially holds that \S\ < \S'\ < 2s. 

It is easy to see that dX C Xy/ U 5', therefore \dX\ < \Xyy \ + |S"| < 6s. The same holds for 
dY. □ 

Applying the above claim we can construct the separator decomposition T for a graph G(V, E) 
as follows. Initially, for the root r of T, let V r = V . For any current node i G T, if 7^ 0, set 
i? = Vj and apply the above claim to get an (X, Y)-separator of S in G[R] with desirable properties. 
Then create two children i and j in T, let Vi = X, Vj = Y, and recursively do the same thing for 
the two children. 

It is easy to see that T is a separator decomposition for G of width at most 6s. For \W\ = 6s, a 
balanced IV"-separator S' can be found in time 2°^ -poly(n) by enumerating all {S, X, Y} partitions 
of W and running the standard network flow algorithm on G[V \ S] to find a separator of X and 
Y. This standard approach for finding balanced separator is also used in construction of tree 
decomposition (see Chap. 11.2 of [29]). Therefore, T can be constructed in time 2°^ -poly(n). □ 

Theorem 14.31 is proved by combining Lemma 14.61 and Theorem 14.51 

5 Counting Algorithms 

This section contains three algorithms for Holant problem with regular symmetric constraint func- 
tions: a simple exponential-time dynamic programming algorithm; a fixed-parameter tractable 
(FPT) algorithm which uses the exponential-time algorithm as a subroutine; an FPTAS on apex- 
minor-free graphs via correlation decay which utilizes the FPT algorithm. 

With the construction of separator decomposition, it is not hard to come up with a very natural 
20(tw(G)) .poly(n)-time dynamic programming algorithm for spin systems by enumerating the vertex 
boundaries and separators of components in the separator decomposition. However, the flexibility 
of Holant problems causes many new issues to the computation, which require more sophisticated 
algorithms to deal with. 

5.1 Simple exp(0(n))-time algorithm 

Any Holant problem can be computed in time exp(0(|£'|)) by enumerating all configurations. For 
Holant problem with regular constraint functions, there is a simple dynamic programming algorithm 
which runs in time exp(0(|V|)). This algorithm is used as a subroutine in our main algorithm. 
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Theorem 5.1. Let = (G(V, E), {f v }vev) be a Holant instance where f v : [g] de sw — > C are 
symmetric functions. If all f v are C -regular for some constant C > 0, then hol(Q) can be computed 
in time (qC)°" v v . 

We enumerate the vertices in V in an arbitrary order v\, V2, ■ ■ ■ , v n . Let G k (V k , E k ) be a sub- 
graph induced by the first k vertices, i.e. Vfc = {vi | 1 < i < k} and E k = {viVj G E | 1 < i, j < k}. 
For a v £ V%, let deg k (v) denote the degree of v in G k - Fix any 1 < k < n. For i = 1, 2, . . . , k, let 
^i*^ be symmetric functions at vertex V{ in the form <$) : [q\ dc &k( v i) — > C. We define the following 
quantity: 



In fact, each Z I k, \<j><£ \ I defines a new Holant problem on G k . And the result of the 

\ I J i=l,2,...,kj 

original Holant problem is given by hol(O) = Z(n, {f Vi }i=i2 ...n)- In general we have the following 
recursion: 



Z(O,0) = 1; 



Z [k, 



i=l,2, 



where 0^ x) 




otherwise. 



This recursion separates the summation into different cases of configurations around v k and mod- 
ifies the functions at the adjacent vertices according to the configuration. The correctness of the 
recursion can be easily verified by observing that the edge set E k is the disjoint union of E^i and 
Ek{v k ). 

We then describe a dynamic programming algorithm which computes the Holant problem 
Z{n, {f Vi }i=i,2,...n) in time (qC)°^ if all f v are C-regular. The algorithm consists of two phases: 

1. Preparation: For every v G V, construct {PlN CT (f v ) \ a G [qY,0 < I < deg(u)} which contains 
all pinning outcomes of f v . For symmetric f v this can be done in time polynomial of deg(u). 

2. Dynamic programming: It is easy to see that for any 1 < k < n and any 1 < i < k, function 

is an outcome of a sequence of pinning of f Vi , Moreover, it holds that 

4f G {PlN CT (/„,) | a G [ g ]deg»(^)-deg fe (^)| > 

where the size of the set is bounded by C since f Vi is C-regular. Therefore, Z I k,\ 4$ \ 

for all 1 < k < n can be stored in an n x C n table, while each q>vi can be retrieved from 
{PlN^ (f Vi ) | a G [q] dc &n( v i)-dcg k (vi) j an j nc j ex ra nging over [C]. It takes at most q°( n ) time 

to fill each entry of the table. The total time complexity is (qC)°( n \ 
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5.2 Fixed-parameter tractable algorithm 

Theorem 5.2. Let = (G(V, E), {f v }vev) be a Holant instance where f v : [gj^sw — > C are 
symmetric functions. If all f v are C -regular for some constant C > 0, then hol(O) can be computed 
in time 2 ( tw ( G )) • poly(n) where tw(G) represents the treewidth of G. 

The 2°( tw ( G )) growth in treewidth is critical to our approximation algorithm on planar graphs 
introduced later, although any faster growth in treewidth is still fixed-parameter tractable. 

The setup. Let U C^bea set of vertices. Let dU denote the vertex boundary of U, i.e. dU = 
{v G V \ U \ 3u G U,uv G E}. Let H(U U dU, F) be the subgraph that F = {uv G E \ u, v G 
[/orfiG[/,ue 917} , i.e. F includes all edges within U and all edges crossing between U and dU 
(but not those edges with both endpoints in dU). For each v G dU, let <j) v : [g] de Sn (") — >. {o, 1} be 
a boolean symmetric function. We define the following quantity: 

z {UAMvedu) = Y II 9v{o-\ F(v) ), (1) 

a£[q] F v€UUdU 

ff v iiveU, 
where q v = < 

\<j) v iivedU. 

In fact Z (U, {4>v} V £9u) defines a Holant problem on graph H(U U dU, F) with function f v at each 
v G U and boolean constraint <p v at each boundary vertex v G dU. And the original Holant problem 
can be written as hol(O) = Z(V, 0). 

The recursion. Suppose that U can be partitioned into S,Ui,U2 such that S is a (C7i, C/2)- 
separator of £7 in G[[7], where G[?7] is the subgraph of G induced by U (note that S is not 
necessarily a separator of H). It is obvious that dU\ C 5U dU and SL^ ^ <S U 9f7, where dU\ 
and 9f/2 are respectively the vertex boundaries of U\ and U2 in G. For each v G £ U <9C7, let 
do{v), di(v),di{v) denote the number neighbors of v in S U c?£7, U±, U2 respectively, i.e. 

do(t;) = \{u G 5 U dU I uv G di(u) = [{« G C7i | uv G d 2 (v) = \{u G C/ 2 | uv G 

It holds that do(v) + d\{v) + ^(v) = deg#(u). 

For any v G 5 U <9f7 and i = 0, 1, 2, each tuple <r G [g]^ 1 ') can be mapped to a boolean function 
Peer ct (g v ) which indicates all tuples that have the same effect of pinning on g v as a, where g v is 
still defined as that g v = f v for v G S and g v = (fr v for v G dU. We call Peer ct (g v ) the peer image 
of <r at u. For each v G 5 U dU and t = 0, 1, 2, let = {Peer ct | cr G [q] di ^} be the range of 
peer images over all a G [gJ^W. Let be a sequence indexed by (j) % v for v G S U d£7 and i = 0, 1, 2, 
such that 0* G V % vl i.e. 0* is a peer image of some a G [gJ^W. 

Due to Lemma [3751 for any w£SU <9£7, for any <jj G [g] d< ("),i = 0, 1,2, the value of g v {p§o\o~2) 
is uniquely determined by the peer images Peer CTj (g v ),i = 0,1,2. For <f> l v G Vl,i = 0,1,2, we 
write g v (<fivi filial) f° r the unique value of g v {o-oo~io~2) for all (<7o> 01,02) with the same peer images 
Peer^ (g v ) = 4 for i = 0, 1, 2. 
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We then have the following recursion for Z (U, {(t>v} v& Qjj) defined in (pQ): 

Z{U,{M v eau)= E Zo(<P) ■ Zi(<P) ■ M<P) U 9*W><f>l<t>l) , ( 2 ) 



<Mie7>i veSUdU 
i=0,l,2 



where Z (cf>) = hoi \H [S U dU] , {$} veSvdu 

z 1 (<f>)=z(u 1 ,{ ( f>l} vedUi ), 
z 2 (<f>) = z(u 2 ,{<f>l} vedU2 ). 

Note that peer images are boolean functions, thus Zq, Z\ and Z 2 are well defined. Also note that 
only Z\ and Z 2 are recursions and Zq is a new well-defined Holant problem which can be directly 
computed by the simple algorithm of Theorem 15.11 

As an example, consider counting matchings, which is a Holant problem of regular constraint 
functions. The peer images cj) l v actually correspond to that vertex v is matched or unmatched^! by the 
corresponding subset of incident edges of v. The Holant problem Zq(4>) = hoi yH [S U dU] , {<fiv} V £sudU 
counts the number of perfect matchings of those vertices that claim to be matched in H [S U dU\. 

We then prove that this recursion holds for Z (U, {4> v } vG q U ) defined in ([1]). 

Proof. Since S is a (U\, [^-separator of U, the edges in the original subgraph H(U U dU,F) can 
be partitioned into five disjoint sets: 



fori = 1,2, Ei 
fori = 1,2, Fi 



internal edges of S U dU, i.e. F = {uv £ F \ u£ S, v £ S U dU}; 
internal edges of Ui, i.e. Ei = {uv G F \ u,v G U}; 
boundary edges of Ui, i.e. F{ = {uv G F | u G Ui,v G 5 U dU}. 



Each vertex v G SUdU has precisely di(v) adjacent edges in Fi for i = 0,1, 2. And for i = 1,2, dUi 
is precisely the set of vertices in S U dU with positive di(v). We can enumerate all configurations 
a G [q] F ° uFlUF2 by enumerating legal local configurations a l v G [q\ Fi ^ for each individual vertex 
v G 5 U dU and each i = 0,1, 2, where being legal means that there exists a a G [q] F ° uFlUF2 such 
that <7 |_f,(j;) = °ti f° r all u G 5 U 9C/ and i = 0, 1, 2. 

For a tuple cr G [q] k , we define the indicator function 1 CT : — > {0, 1} as that 1 ct (t) = 1 if and 
only if t = cr. Then we can rewrite ([I]) as follows: 

zZ ^(^'{^I^wj^ac/J '^(^'{^wl^J II S« (ff k(«)u^(«)uJ!,(«)) 
= ^ ^(^.{VJ^)-^^,^}^) II 9v {°VvVv) ■ (3) 

legal <Tje[g] F iM i>£SU<9C/ 
i=0,l,2 

In fact, Fi can be partitioned into disjoint -Fi(v) for v G 3£/i and i*2 can be partitioned into 
disjoint ^(f) for v G c5C/2- Thus for i = 1,2 all local configurations {cr* G [fj^^j-vesudc/ are legal. 



3 We can ignore the overmatched cases in our discussion because they nullify the configuration. 
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A collection { a v} v ^su8U °^ ^ oca ^ configurations of edges in Fq of individual vertices is legal if and 

has value 1 (it has only two possible 



v£SUdU 



only if the Holant problem hoi [H[S U dU], {l a o } 
values or 1 as every indicator function has value 1 on exactly one input). Thus we have 



Wv€SUdU 
i=0,l,2 



v&dU 2 



J] 9v(°°M)- (4) 

vesudU 



For v e S U dU and i = 0, 1,2, fix 0* G "P*, i.e. 0* = Peer ct for some a £ [q] di ^ v \ We can 
group configurations {c l v } v£ sudu i=o l 2 into equivalence classes {cr* € [al^W | Peer ; (g v ) = (j) 1 } 
according to their peer images. Due to Lemma 13. 3[ configurations from the same class yields the 
same value of g v (a®al<7%) = g v (4>v, <f>l, 4>D- Therefore, 



®= E 



E 



a°&{ q ] F a(v):P W R(* v ,g v )=<P v 
\/v£SUdU \ Vv^SUdU 
i=0,l,2 



/ 

V 



E 



o-j6[g] i? iW : PEBR( < Ti,g„)= 

Vvesudu 



E 



^ v elq] F 2M:PEER(alg v )-- 
VveSUdU 



hol[H[SUdU],{l a o} veSudu 
J 



\ 



I 



\ 



J 



n u<t>i<t>iAD- (5) 

vesuau 



And any peer image (j> l v € V\ is a boolean function which indicates all such a G that have 

the same peer image Peer ct (g v ) = <j) l v , thus it is straightforward to verify the following identities: 



h6L{H[S\ja^,{^} vegugu 

z ( u id^Uau 1 
z ( u ^{^Uau 2 



E 


hol(H[SUdU],{l a o} veSudu ) 


1 '):Peer( ( tS, 9 „) = 

VveSuau 


i>°v 


E 


Z ( U ^{ 1 °l} v eau 1 )> 


■"):PEER(o-J,g t ,)= 
VveSltdU 


Pi 


E 


Z { U ^{^ v } veaU2 ). 


u >:PEER(o*,ff„)= 
WvGSUdU 


Pi 



Substituting these identities back in ((5|), we deduce the recursion ([2]). 



□ 



The algorithm. We then describe an algorithm which computes hol(G(V, E), {f v }vev) hi time 
2 0(tw(G)) . po ly(|y|) if all f v are C-regular for some constant C > 0. 



18 



1. Constructing separator decomposition: By Theorem 14.31 a separator decomposition T of 
input graph G of width 3(tw(G) + 1) can be constructed in time 2°( tw ( G )) • poly(n). 

2. Enumeration of peering closures: For every v £ V and each < k < deg(v), construct set 
P(v,k) = {Peer ct (/„) | a £ [q] k } and all possible unions (defined on boolean functions) of 
members of P(v,k). Due to Lemma 13.61 this will cover all functions in Peer* (f v ). And 
due to Lemma 13.71 there are at most 2 C possible unions. The total time cost is polynomial 
because for symmetric functions all such P(y, k) can be constructed in polynomial time. 

3. Dynamic programming: Let T be the separator decomposition constructed in the first step. 
Then each node i £ T associated with a vertex set Vi and a separator Si such that \Si\ = 
0(tw(G)) and \dVi\ = 0(tw(G)), and if j and k are the two children of i in T, Si is a (Vj, V k )- 
separator in G[V^]. Apply the recursion (J2]) in this tree structure as follows: For each leaf 
I £ T, Vi = 0, and Z(Ve, 0) = 1; and for each non-leaf node i £ T with children j and k in T, 
Z(Vi, {4> v } V £dVi) is computed according to the recursion ([2]) by setting U = Vi, U± = Vj and 
U 2 = V k ; in particular for the root r of T, V r = V and Z(V, 0) = hol(G(V, E), {f v }vev)- 

There are 0(|V|) nodes in a separator decomposition. For all Z(Vi, {<j)v}v&dVi)i every <j> v 
is a boolean function in Peer* (f v ). Due to Lemma 13.71 since /„ is C-regular, once Vi is 
fixed there are at most 2 C possible <f> v for each v £ dVi, where \dVi\ = tw(G). Therefore all 
Z(Vi, {4>v}vedVi) can be stored in a 0{n) x 2°( c '- tw ( G )) table. 

Each entry of the dynamic programming table is filled according to the recursion ([2]) , which 
involves three nontrivial tasks: 

(a) (computing Zq): Due to Lemma 13.71 an Y 'Pv ^ Peer* (f v ) is still C-regular since /„ is 
C-regular, thus Zq = Yvo\{Hi\Si U3Vi], {4>v}v€Siii8V^) 1S a Holant problem with C-regular 
constraint functions which can be computed in time (qC)\ SiUdVi \ = (qC) c, ( tw ( G ')) due to 
Theorem 15.11 

(b) (evaluating g v ): Each g v (c^, (frl, 4>1) can be easily evaluated by evaluating g v {(jQ(J\cr2) for 
arbitrary £ <p® , o\ £ <fil , a 2 € 4>1 . 

(c) (computing the sum): For every v £ Si U dVi, enumerate all < 2 C possible boolean 
functions of appropriate arity (f> v £ Peer* (f v ). The total time is bounded by 2°( c ' tw ( G ^ 
because |5» U dV\ = 0(tw(G)). 

The time cost for filling one entry of the dynamic programming table is bounded by 2°( tw ( G )) 
for constant C and q. 

The total time cost for the above algorithm is bounded by 2°( tw ( G )) • poly(re) for constant C and q. 
5.3 FPTAS from correlation decay 

Theorem 5.3. Assume the tractable search for the Holant problem Holant(<5, J 7 ) where Q is an 
apex-minor-free graph family and T is a regular family of nonnegative symmetric functions. The 
strong spatial mixing implies the existence of FPTAS for Holant(£/, J-). 

Let f2 = (G(V,E), {f v }vev) be a Holant instance, where G is an apex-minor-free graph and all 
fv '■ [q] dcg ^ — > M + are C-regular symmetric functions for some constant C > 0. Let [i be the Gibbs 
measure defined by the Holant instance £1. 
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Assume the tractable search and strong spatial mixing for Holant (£/, F). We have the following 
lemma for approximation of marginal probabilities. 

Lemma 5.4. Let e G E. Let ACE and t\ 6 [q]^ be a feasible configuration. The marginal 
probability nl A (i) for any i G [q] can be approximated within any additive error e in time poly(n, -). 

Proof. Let N r {e) = {e! G E | dist(e, e') < r} be the r -neighborhood of edge e in G. Let B r (e) = 
{uv G E \ N r (e) | 3wv G N r (e)} be the vertex boundary of the r-neighborhood. 

Denote A = B r (e) \ A. When the strong spatial mixing holds, by Definition 12.41 f° r an Y 
<ta,7ta G [q] A that both (ta,<ta) and (ta,7ta) are feasible, it holds that \\[il A,rJA — \i e K ^^ || TV < 
poly(n) • exp(f2(— r)). Therefore for any <ta G [q] A that (ta,o"a) is feasible, we have 

II^ta _ ^a^a|| tv < poly(n) . exp(n(-r)), (6) 

because is a linear combination of all such [i A '° A . 

Note that the joint configuration (ta,o"a) fixes the boundary B T (e). Thus for each i G [q], the 
marginal probability fil A ' aA (i) can be computed precisely from the r-neighborhood as follows: 

Let W be the set of incident vertices of N r (e) and F = N r (e)\A. Let H(W, F) be the subgraph 
formed by removing edges fixed by a\ from the r-neighborhood. For i G [g], let e i — ?- z denote the 
configuration on {e} that simply assigns value i to edge e. We have that 

hol(H'(W,F\{e}),{f^ e ^\ ) 

^ e lj hol(H(W,F),{fv A ' aA } veW ) ' [) 

where f T v = Pin Tij (f v ) and r v = r . 

The correctness of the equation and the well-defined-ness of the new Holant problems are easy to 
verify. 

Since the original graph G is apex-minor-free, due to Theorem 12.61 we have tw(H) = 0(r). 
Since all original /„ are C-regular, then trivially all f% are C-regular since they are just results of 
pinning f v . Then applying Theorem 15.21 the new Holant problems defined in ([7]) can be computed 
in time 2°^ ■ poly(n). Therefore, the marginal probability with boundary condition He A,tTA (i) for 
any i G [q] can be computed precisely in time 2°^ ■ poly(n) once a feasible (ta,o"a) is given. 

Due to the tractable search for Holant (G, J 7 ), given any feasible t\ G [g] A it is possible to 
efficiently choose an arbitrary feasible a G [q] E agreeing with t\. Thus a feasible <ta G [q]^ can be 
efficiently constructed by restricting the aforementioned a on A. Due to Q, the original marginal 
probability nl A (i) for any i G [q] can be approximated within an additive error e in time poly(n, -) 
by choosing appropriate r = 0(logn + log -). □ 

With the above lemma, we can apply the standard self-reduction procedure to obtain the FPTAS 
for Holant (Q, F). 

Let r G [q] E be a feasible configuration, i.e. the Gibbs measure /i(r) > 0. Enumerate edges in 
E as ei, e2, . . . , e m . For each < k < m, let Ej, = {e±, . . . , e^}, G [q] Ek be consistent with r on 
Ek, and pk = ^ _1 (r(efc)). The following identity hold for h(t): 

m m 

fi(r) = ] Pt e [a(e k ) = r(e k ) \ cr(ei) = r(ej), 1 < i < k - 1] = JJ p k - 



k=i 
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On the other hand, /x(r) = Uvev ^\ B{v ^ ■ Thus hol(O) = ^fij^f • If for each ( X ) P* can 
be approximated in an additive error e; and (2) p^ > is a constant, then the product n^LiPfc can 
be approximated within a multiplicative factor (1 ± 0(ne)). While (1) is guaranteed by Lemma 
15.41 (2) can be achieved by trying /ie^ -1 ^) for all i G [q] and choosing r(efc) to be the i with the 
largest returned value. This gives us an FPTAS for the Holant problem. 

We then can directly apply any known strong spatial mixing result to get the FPTAS. For 
example, combining with the result of [35], we have the following corollary. 

Corollary 5.5. There exists an FPTAS for counting q-coloring on apex-minor-free triangle-free 
graphs of maximum degree at most A if q > aA — 7 where a ~ 1.76322 is the solution to a a = e 
and 7 = 4« 3 -W-3a+4 _ 0.47031. 

Note that although the original result of |35j is proved for single-site strong spatial mixing where 
the boundaries differ on only one vertex, it implies our definition of strong spatial mixing on any 
finite graphs. 



6 Correlation Decay 

In this section we apply the recursive coupling technique |35| to Holant problems, and prove strong 
spatial mixing for subgraphs world |36| and ferromagnetic Potts model. The algorithmic implica- 
tions of these correlation decay results are presented in the end of this section. 



6.1 Recursive coupling on Holant Problems 

Consider a Holant problem Holant(£?, J-) and an instance f2 = (G(V, E), {f v }v£v)- Let R C E, 
called region. Define 5R = {uv G E \ uv R, 3uw G R} the edge boundary of R. Define 
V R = {v G V I 3uv 6 R}. 

A boundary configuration of R, is a a G [q] SR - For every boundary configuration a G [q] SR and 
a configuration 77 G [q] R of the region R, define the regional weight as 

w r(v) = ~[l fv(v \r(v) tr \(6R)(v)), 
veVn 

where 77 is the restriction of r/ on the edges in R incident to v, a |(<5_r)(„) is the restriction of a 
on edges in 5R incident to v, and f v (a \jir v \ r\ \(sr)( v )) evaluates f v on the concatenation of them. 

We say that a boundary configuration a G [q] SR is R-feasible if there exists an n G [q] R such 
that w R (rf) > 0. For i?-feasible boundary configuration a G [q] SR , a regional Gibbs measure \i R over 
[q] R can be defined as that fJ. R {r]) = — Wr ^ v \,\ for each n G [q] R . For R' C R, let fJ. RR , denote 
the marginal distribution of \x a R on R', and we write that \x a R e = fi R ^ e y 

Definition 6.1. Let R C E, R' C R } and ct,t G [q] SR be two R-feasible boundary configurations. 
Let *&(R, o~, r) be a coupling of fi R ,fi R . Define the discrepancy of ^(i?, <r, r) on region R' C R as 

Disc* (Ri(T)r ) (Of) = Pr [rj \ R >^ 7} \ R ,] 
If R' = {e}, we write that Disc^ (iW) (e) = Disc^ {jRj(Tir ) ({e}). 
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Definition 6.2. For any two R-feasible boundary configurations o~,t G [q] SR differing on A C <5i£, 
a sequence of R-feasible boundary configurations a\, o~2, ■ ■ ■ , o~t is called a feasible path from a to r 
if a = (Ti,t = at and o"j,crj +1 differ only at one edge e G A for each 1 < i < t. Let T(a,r) be the 
minimum such t, or be oo if no such path exists. 

Lemma 6.3. Let A C E and cr, r G [g] A be two feasible configurations differing on A C A. Lei 
R = E \ A and e G -R. There exist two R-feasible boundary configurations o\t' G [g] 5R differing 
only on edges in A smc/i i/iai 



I Me MellTV 



< T{a\ t') ■ max ||^ - /^Jtv 

TV (Tl,(T2E[?] ia 

differ on an e'sA 



< T(cr / , t') ■ max Disc w(Jij(J1)(72 ) (e) , 

differ on an e'GA 

/or arbitrary coupling ^(R, a\, a 2) of /i^ 1 , /i^ 2 . 

Proof. Let cr', r' G [q] SR consistent with a, t on 5R respectively. It is easy to check that a', r' satisfy 
the equation. Let a' = o~\, o~2, ■ ■ ■ > o~t = r' be the feasible path from a' to r' of length t = T(a' , r') 
and \i a Re — [i T Re can be bounded by applying path coupling to /J^ e , /J^ 1 for 1 < i < T{a' , r'). 
The last inequality is due to the coupling lemma. □ 

The above lemma reduce the strong spatial mixing to the discrepancy witnessed by a coupling 
of /i R , jj, T R with o", r disagreeing at one edge. We then show a way to recursively construct the 
coupling. This method is proposed by Goldberg et al. in [35] on colorings. 

The recursive coupling. Let R C E and eo G R. Let a, r G [a]' 5 ^ be any two i?-feasible 
boundary configurations that differ at only one edge e G SR. Let -R(e) be set of edges in R 
incident to e. Let ^R( e ) (R, o~, r) be a coupling of marginal distributions /i R R ^ , \i r R R ^ . A coupling 
^f(R, cr, r) of regional Gibbs measures /j r , /i R can be recursively constructed by the local coupling 
rule ^R( e }(R, a, r): Let (??,?/) G [g] fl x [g]^ denote the pair sampled from ^(R, a, r). 

1. (Base case) If eo G -R(e), sample (77 |i?( e ),7/ |_R( e )) according to ^We) (-R, c, r) and arbitrarily 
sample the rest of (77, 7/) conditioning on (7/ |fl( e ), 7/ U(e)) as l° n g as (?7j ?/) is a faithful coupling 
of /j, r , n R . If |-R(e)| = 0, in which case e and eo are disconnected in G[Vr], sample (77,7/) such 
that (r/(eo), V( e o)) is perfectly coupled. 

2. (General case) |i?(e)| > and eo G" R(e). Sample (77 |_R( e ),r/ U( e )) according to ^ji( e \(R, a, r). 
Let x,y G [g] R ( e ) be two configurations that x = rj |_R( e ),y = 7/ U(e)- Construct new region 
and boundaries as: i?' = R\R(e); a x G [g]* 5 ^ agrees with <r on common edges and a x |_R( e )= x; 
and r' y G [q] SR ' agrees with r on common edges and t v |_R( e )= V- The rest of (r],rf) is sampled 
from a coupling ^>(x, y) of /i^, , \i Rl . If x = y, then = t v and ^(x, y) is a perfect coupling. 
If x 7^ y, let cr^ , . . . , 0-£^ be a feasible path from a x to r y of length t = T(a x ,r y ). Let 

y) be the composition of coupling \I> ^i?', cr^'^, cr^'^J , i = 1, 2, . . . , t — 1, in the same 

manner as path coupling, where each ^ ( i?', cr^'^ , cr£? ^ ) can be recursively defined as a'f^ 



'i+l 



and differ at only one edge. It is easy to verify that ^(x,y) is a coupling of fi R , n R 

This complete the construction of ^{R, a, r). 
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The following lemma is similar to the one proved in [35] for the recursive coupling constructed 
on spin systems. 

Lemma 6.4. For the coupling ^>(R, a, r) constructed as above, we have 
Disc^ (iW) (e ) 

T(a x ,ry)-1 

xjty 

< Disc* KCe)(jRi(T , r ) (e) • max : T^, <t 2 ) • m ^ Disc^ (RVR(e)j(TliCT2) (e ) 

differ on R(e) differ on an edR(e) 

Proof. The lemma follows directly from our construction of ^(R,a,r). □ 

A standard choice of ^_R( e ) (R, o~, r) is the one that Pr^y)^^,,^^,-) [77 \r( b ) = v' \r(b)] 1S max- 
imized, i.e. Disc<j> (r,ct,t) (R( e )) i s minimized. 



6.2 The subgraphs world problem 

The subgraphs world model used in [36] for developing FPRAS for the ferromagnetic Ising model, 
is a counting problem computationally equivalent to the Ising model under holographic transfor- 
mation. 

Definition 6.5 (Subgraphs world [36]). The subgraphs world with parameters (A, jj) defined as 
follows. Let G = (V, E) be an undirected graph. The subgraphs world partition function is defined 
as: 

Z suh (G) = £ ^WIaW, 

XCE 

where odd(X) denotes the set of vertices with odd degree in the subgraph (V, X) . 

The subgraphs world with parameter (A, jj) can be interpreted as a Holant problem on incident 
graph Tq as follows. The incident graph Tq has left vertex set V and right vertex set E, and for each 
v € V and e G E, (v, e) is an edge in Tq if e is incident to v in G. The function on each left vertex 
v is [1, /i, 1, fi, . . .] and the function on each right vertex e is [1,0, A]. Let be the Holant instance 
defined as above. It is easy to verify that all functions in £1 are 3-regular and Z SU \>(G) = hol(O). 

For convenience of analysis, we consider the following equivalent Holant problem which is defined 
on the original graph G instead of the incidence graph. Let £1' = (G(V,E), {f v }vev) be a Holant 
instance where each f v = [/o, /i, . . . , fdeg(v)] has that fk = ii\ k ^ 2 if k is odd and fk = A fc / 2 if k is 
even. It is easy to see that hol(O) = hol(f2') and also the two Holant problems have exact the same 
Gibbs measure. Thus it is sufficient to analyze the correlation decay on this new Holant problem. 

Theorem 6.6. Let Holant (Q, J 7 ) be defined by the subgraphs world of parameter (ft, A) with < 



/j, A < 1 on graphs with degree bound A. If A < ^jz^r~> then Holant {Q, F) has strong spatial 
mixing 
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Proof. Let R C E be a region and <r, r € {0, 1} SR be two i?-feasible boundary configurations 
differing at e G 5i? satisfying cr(e) = and r(e) = 1. The regional weights w R ,w R and regional 
Gibbs measures //^ can be defined accordingly. 

Let ^ R(e) {R,a,T) be the coupling of ^ R{e) ,^% R{e) such that Pr M) ^ R{e){R ^ T) [r] \ R{e ) = 
r( t^(e)] is maximized, i.e. the discrepancy Disc^, H{e) (/j j(Jjr ) (R(e)) is minimized. We first give an 
upper bound on Disc# ( Bj0 . iT ) (-R(e)). 

Let i] € {0, 1}^ be a configuration on R, we use to denote the number of edges in R(e) that 
are assigned to 1 by ij. Let 



^€{0,1}^ r;e{o,i} fl 

n^is even n v is odd 

Then for every r/ E {0, 1}^, we have /^(t?) = - "j_ T ' w and 

— " ',' 11 i , if n„ is even 

W e -fJ,+W /fl ' v 

It holds that 

■r ff / n r/ N ^r, «V /J W (l - p 2 ) 

if is even, fi R (r]) - ^ R {rj) = ■ ■ — = w v ■ ■ — > 0; 

W e + W W e - [i + Wo/H (W e + W )[W e ■ fl z + W ) 

if n v is odd, - fifa) = = ■ < + W if~ l \+ \ < °" 

«; e + w; w e - fi + w /fi [w e + w )(w e ■ ^ + w Q ) 



Thus in the coupling V^Rfe) (R, &, t) , we have 



u> (l - /U 2 ) 



»^ e) (R,«,r)(R(e))= Ml) ~ Ml) = E ^■ K + Wo)(We . M 2 + Wo) 

■q: n rj is even 7/: rirj is even 

w e tt; (l - fi 2 ) 1 - /j. 2 



(lO e + W ){w e ■ fl 2 + W ) (1 + W /w e )(l + /J 2 • W e / Wo) 



(8) 



We then show that A < ^ < ^-y. We assume a total order on all edges. For any 77 with even n v , 
let ^(77) be the configuration resulting from flipping the state of r\ on the first edge in R(e). Note 
that 4> is a bijection between configurations in {0, 1}^ with even n v and those with odd n v . It is 
easy to check that 

if n v is even, w v > Xw^y, and if n v is odd, w v > A/j 2 • w^-it^y 

Combining with the fact that is a bijection, we prove that A < ^ < ^j. Substituting this 

into (jSJ), we have Disc^ fi(e) (/j iCTT ) (R(e)) < r^^rp ■ And since /x, A > 0, all boundary configurations 
are i?-feasible, thus for any boundary configurations a', t' differing on t edges, we have T{a' , r') = t, 

i.e. we can migrate from one boundary configuration to another by modifying one edge at a time 

( 1+Xu 2 ) 2 

without violating the feasibility during the process. Therefore, if A < v -,_ 2 , then 

Disc^ (e)(B)CT)T) {R(e)) ■ ^T^^Tfaa*) < ^~^ 2)2 ■ \R(e)\ < 1 • A = 1. 

differ on -R(e) 
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For any eo £ R, we can apply the recursion in Lemma [63] for dist(e, eo) many times where dist(e, eo) 
denotes the distance between e and eo, thus Discern atT \ (eo) = exp(J7(— dist(e, eo))). Then applying 
Lemma |6, 31 since T(a',r') = poly(n) for any ct',t', we have the strong spatial mixing. □ 

The ferromagnetic Ising model is a spin system specified by <&e '■ {0, l} 2 — > K + and $y : 
{0, 1} — > IR + such that &e(,x, y) = a> \\ix = y and <&e{x->v) = 1 if other wise;, and $y(x) = b > 
if x = 1 and $v(x) = 1 if x = 0. We call (a, 6) the parameters of the system. 

The Ising model can also be specified by the inverse temperature (5 and external field B as 
follows. Given a graph G = (V,E), the partition function is defined as 

Zj sing (G) = Y, II exp(/3a(z)a(j)) J] exp(Sa(z)) 
<7e{-i,i} v ' *je£ iGV 

Theorem 6.7 (Jerrum and Sinclair [Mj)- G(V, E 1 ) 6e a graph. Let X = = tanh/3,/u = 
Ijqrj-I = tanhi?, then 

Rising (G) = Mg • Z suh (G), 
for some Mq which can be computed in polynomial time. 

The transformation from the ferromagnetic Ising model to the subgraphs world model is actually 
a holographic transformation and the above theorem can be seen as a special case of Valiant's Holant 
theorem [131156] . 

Translating the conditions in Theorem 16.61 for subgraphs world back to the ferromagnetic Ising 
model, we have that A < ,) y2(h or equivalently A < \ 2B , 2 . , 2 , 2U ' . 



6.3 Recursive coupling on spin systems 

Let G(V,-E) be an undirected graph, and $ : [g] 2 — > M + be a symmetric function of nonnegative 
values. Consider the g-state spin system whose partition function is defined by 

Z(G) = ^2 w ( a ) where w(a) = $(cr(u),a(v)). 

vG[q] v {u,v}eE 

With this definition of weight, for any feasible a £ [q] for A C V, we can accordingly define the 
Gibbs measure fi a over [q] v and the marginal distribution at jj% at vertex v. 

Definition 6.8 (Strong Spatial Mixing on Spin Systems). A spin system on a family of graphs has 
strong spatial mixing (SSM) if for any graph G = (V,E) in the family, any v £ V, A C V and any 
two feasible configurations o~\ 1 t\ £ [g] A ; 

\\K, A -M; a Htv < exp(-ft(dist(t>,A))), 

where A C A is the subset on which a\ and t\ differ, and dist(i> , S) is the shortest distance from 
v to any vertex in A. 

It is easy to verify that under this definition SSM of a spin system is equivalent to the SSM of 
the Holant problem on the bipartite incident graph Xq which simulates the original spin system. 
Thus to use the FPTAS which rely on the SSM for Holant problems, it is sufficient to prove SSM 
on the original spin system. 
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Let R C V, called a region. Let dR = {v G V \ R \ 3uv G E,u G R} be the vertex boundary of 
R and 5R = {uv G E \ u G R,v i?} be the edge boundary of -R. Denote by E(R, R) = {uv G E | 
u £ R,v £ R} the set of internal edges in R. 

An edge boundary configuration of R, or just boundary configuration for short, is a ex € [q] SR - 
For every boundary configuration a G [g]^ and a configuration rj G [g;]^ of the region i?, define the 
regional weight as 

w%(??)= n(v)) ■ \{ *(o-(e), tj(u)). 

uv£E(R,R) e=uv£SR 

uGR 

We say that an edge boundary configuration <r G [g]* 5 ^ is R-feasible if there exists an 77 G [g] K 
such that w R (r]) > 0. For ^-feasible boundary configuration a G [9]^, a regional Gibbs measure 

UK over [g] R can be defined as that [i R (r)) = ^ — "^lu 17 (tt) ^ or eacn ^ £ feP- For v £ R, let be 
the marginal distribution of /i^ at vertex v. 

Definition 6.9. Let R C V , v £ R, and a, r G [q] SR be two R-feasible boundary configurations. 
Let ^{R,a,r) be a coupling of fi R , fi R . Define the discrepancy of ^f(R,a,r) at vertex v as 

Disc«r (Ji)CT)r ) (v) = Pr [v(v)^ r)'(v)] . 

Lemma 6.10. Let Ac^ and 0", r G [g] A be two feasible vertex configurations differing on A C A. 
Let R = V \A, v £ R and V = {uv G J5 | v G i?, u G A} &e i/ie set 0/ boundary edges of R incident 
to vertices in A. There exist two R-feasible edge boundary configurations o~',t' G [q] SR differing 
only on edges in V such that 



IItv 



< T(ct / , r') • max " ^r,Jtv 

differ on an eSV 



< T(a', t') ■ max Disc^ {jRj(Tli(72) , 

o"l,o- 2 G[(?] dK 
differ on an e€V 

/or arbitrary coupling ^(i?, o"i, 02) of fi^ 1 , fi R . 

Proof. Let a', t' G [g] 5 ^ be defined that the assigned value of each boundary edge is consistent with 
that of the incident boundary vertex in A assigned by a, r respectively. The rest is the same as 
proof of Lemma 16.31 □ 

We describe the recursive coupling for spin systems introduced in |35| . 

The recursive coupling. Let R C V and vq G R. Let a, r G [q] SR be any two i?-feasible 
boundary configurations that differ on only one edge uv G 5R where v G R and u R. Let 
Er(v) = {wv G E I w G R} be set of internal edges of R incident to v. Let ^ v (R,a,r) be 
a coupling of marginal distributions /^j,,, at vertex u. A coupling ty(R,a,r) of fi R , fi R can 
be recursively constructed by the coupling rule *& V (R, a, r) at the vertices f incident to the only 
disagreeing edge: 

Let (r),rf) G [q] R x [g] R denote the pair of configurations sampled from ^(R, a, r). 
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1. (Base case) If v = vq, sample {n{v ), n'(v)) according to ^f v (R, a, r) and arbitrarily sample the 
rest of (77,7/) conditioning on (n(v ), n'(v)) as long as (n, n') is a faithful coupling of fi R , /j, r . If 
v ^ vq and \Er(v)\ = 0, in which case v and vq are disconnected in G[R], sample (t?, t/) such 
that (r](vo),r}' (vq)) is perfectly coupled. 

2. (General case) v 7^ vq and v has neighbors in R, i.e. |-Er(i;)| > 0. Sample (n(v ), n'(v)) 
according to ^ V (R, a, r). Denote (x,y) = (r](v), rf (v)). Construct new region and boundaries 
as: R 1 = R\ {v}; a x G [q] SR agrees with a on common edges and assigns x to edges in Er(v); 
and r x y G [(/] <5 " R agrees with r on common edges and assigns y to edges in Er(v). The rest of 
(rj, rj') is sampled from a coupling ^(x, y) of \i R i , /j^, . If x = y, a x = t v and ^(x, y) is a perfect 
coupling. If x 7^ y, let er^ , . . . , er^'^ be a feasible path from to r 3 ' of length t = T (a x ,T y ). 
Let ^(x, y) be the composition of coupling ^ ^i?', a^ x ' y \ cr^^^j , i = 1, 2, . . . , t— 1, in the same 

manner as path coupling, where each ^ ( i?', of* , ) is recursively defined. It is easy to 



verify that ^(x,y) is a coupling of fi R ,,i/ Rl . This complete the construction of ^(R, a, r). 
The following lemma is proved in [35] for the recursive coupling constructed as above. 
Lemma 6.11 (Goldberg-Martin-Paterson [55]). 

T(o- x ,T y )-l 

Disc* (JW) («o) < 2 Pr hW^A^)^]- £ Disc *(fl\{,}, CT f-),^f)) ^ 

< Disc*,(i? )CT)T ) (v) • max T (a u a 2 ) • ma * DisCv &(fl \ {tl}i(71i(T2 ) (u 

differ on Er(v) differ on an e£En(v) 

6.4 Ferromagnetic Potts Model 

The Potts model is that a g-state spin system defined as that <I>(x, y) = A if x = y and <l>(x, y) = 1 if 
otherwise. We also write A = e' 3 where /3 is the inverse temperature. The weight of a configuration 
<7 G is w{o~) = \ mon <?( E ) 5 where where mon ff (£) is the number of monochromatic edges, i.e. edges 
uv G E that (j{u) = cr(v). A Potts model is ferromagnetic if (3 > or equivalently if A > 1. 

Theorem 6.12. Let Q be any family of graphs whose maximum degree is bounded by A. A ferro- 
magnetic Potts model on graph family Q has strong spatial mixing if q — 2 > (A — 1) (A — 1) A A , or 



in terms of inverse temperature if (3 < 



A+l 



Proof. Let R C V be a region and a,r £ [q] SR be two -R- feasible boundary configurations differing 
on e = uv G <5i? with v £ R and u ^ R. The regional weights u^, and regional Gibbs measures 
A 4 ??) A 4 ]? can be defined accordingly. 

Let ^ v (R,o-,t) be the coupling of fi Rv , fi Rv that Pr^ ,,/)~<i-„(_r i( 7,t) [v( v ) = 7 /( t 0] i s maximized, 
i.e. the discrepancy Disc^(^ i0 - iT ) (v ) at t> is minimized. We first give an upper bound on Discj,^ >rT>T \ ( 

Consider a boundary configuration a 1 that agrees with a on 5R\{e}. Let cr'(e) = g, a free color 
not in [q], and override the definition of function <I> such that $(i,q) = &(q,i) = 1 for all i G [q]. 
We let Q = X^eUp w r (v)- Without loss of generality, assume <r(e) = 0, r(e) = 1 and Co > c\. 

rj(v)=i 
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Denote that c = Yli=2 c * ■ We nave 

cq+Aci+c ' 



j — n. _ f Aci ■ _ -i . 



I Aco+ Cl+C ' 1S»<9, [ C0 +Ac 1+ c' *-U,2<l<g. 

r;(i))=? r](v)=i 

For ferromagnetic Potts model, A > 1 and Aco + ci + c > Co + Aci + c, thus J2rie[q] R ' u ??( 7 ?) < 
J2rie[q] R ^i?( 7 ?) ^ or au ^ 7^ 0. Therefore in the coupling ^ v (R,a,r), we have 

rj(v)=i 



Disc^R^r) (v) = E /^fa) ~ E 



Ac c 



Ac + ci + c c + Aci + c 

r/£{q] R V e[q] R 
r](v)=0 7](v)=0 

c (A-l)(ci +Aci + c) < (A-l)c 



(Ac + ci + c)(c + Aci + c) c + Ac + c 



The value of the last term J^^^ depends on R, a and e. For fixed e, define v(R, a) = C( ^ A *p+ c i 
where Co, c are defined by R, a, and e as above. 

Let R' C i? be any region containing u . Let 7r € [q] SR be an iZ'-feasible boundary configuration 
that agrees with a on common edges and maximizes v(R',ir). The next lemma is very similar to 
the one proved in |34j . 



Claim 6.13. v(R,o) < v{R',ir). 

Proof. Let rj be an .R'-feasible boundary configuration of R' such that 77(e) = q the free color and 77 
agrees with a on all other common edges. We use cj to denote the sum of weight of configurations 
in R' which assign spin i to v with boundary configuration 77 and c v = Yll=2 c i- R- ecau that 
v(R,a) = \ + i+l/ CQ ■ Then the claim follows from the fact that c/co is a convex combination of 
c^/cq over all such rj. □ 

Consider R' = {v}, suppose v has k neighbours in R U dR, where < k < A {u is not taken 
into account as it is always assigned free color in the definition of qs). We want to find a boundary 
condition tt that maximize u(R',tt), which is equivalent to minimize c/cq. 

Since A > 1, the ratio achieves its minimum when tt assigns all the k incident edges of v with 
spin 0. In this case, Co = A fc , c = Yll=2 c i = Q ~ Thus assuming that q — 2 > (A — 1) (A — 1) A A , 
we have 

f A — l)A fc (X — 1)X A 1 

M-M^l (») < v(R,a) < v(R',n) < - _\ - ( > + < - _\ - ^ ^ < j. 

Note that for the ferromagnetic Potts model, all configurations are feasible, and \Er(v)\ < A, thus 

ma £,w T (01,02) = \Er{v)\ < A. 

differ on Er(v) 

Applying Lemma 16.111 and Lemma I6.10| we have the strong spatial mixing. 

□ 
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6.5 Algorithmic implications 

Both subgraphs world and ferromagnetic Potts model are Holant problems of regular constraint 
functions and both satisfy tractable search. Then by Theorem l5.3l we have the following algorithmic 
results. 

Theorem 6.14. Let Q be the family of apex-minor-free graphs of maximum degree A. 

• If A < ^-fz^r-j there exists an FPTAS for subgraphs world of parameters < fj,, A < 1 on 
graphs from Q. 

. If A < 

g2s ( e 2/5 + ^ 2 (^2s + 1 j 2 , there exists an FPTAS for ferromagnetic Ising model of inverse 
temperature /3 and external filed B on graphs from Q. 

2 / q— 2 \ 

• If j3 < — > there exists an FPTAS for q-state ferromagnetic Potts model of inverse 
temperature (3 on graphs from Q. 

The FPTAS for Ising model is not by applying Theorem 15.31 but due to the transformation 
between subgraphs world and Ising model. 
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